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We present the ground and excited state spectra of doubly charmed baryons from lattice QCD 
with dynamical quark fields. Calculations are performed on anisotropic lattices of size 16® x 128, with 
inverse spacing in temporal direction = 5.67(4) GeV and with a pion mass of about 390 MeV. A 
large set of baryonic operators that respect the symmetries of the lattice yet which retain a memory 
of their continuum analogues are used. These operators transform as irreducible representations 
of SU(3 )f symmetry for flavor, SU(4) symmetry for Dirac spins of quarks and 0(3) for spatial 
symmetry. The distillation method is utilized to generate baryon correlation functions which are 
analysed using the variational fitting method to extract excited states. The lattice spectra obtained 
have baryonic states with well-defined total spins up to | and the pattern of low lying states does 
not support the diquark picture for doubly charmed baryons. On the contrary the calculated spectra 
are remarkably similar to the expectations from models with an SU(6)x 0(3) symmetry. Various 
spin dependent energy splittings between the extracted states are also evaluated. 

PACS numbers: 12.38.Gc, 14.20.Mr 


I. INTRODUCTION 

The study of hadrons containing charm quarks has recently undergone a renaissance. This resurgence of interest 
started with the discovery of new resonances in the charmonium system as well as a few charmed baryons. So far, 
emphasis has been given to the study of the meson sectors both theoretically and experimentally, while heavy baryon 
physics has received substantially less attention. Similar to charm mesons, a comprehensive study of charm baryons 
can provide similar insight into the strong interaction. However, in comparison to the many light and strange baryon 
states, only a handful of charm baryons have been discovered and a reliable determination of the quantum numbers 
of most of these observed charmed baryons has yet to be made [T]. Only very recently a few excited singly charmed 
baryons were discovered. Similarly, there is no observation for triply-charmed baryons although QCD clearly predicts 
such states. For doubly-charmed baryons, only SELEX reported the discovery of five resonances and interpreted those 
as S+^(3443), 5+^(3520), S++(3460), S++(3541) and S++(3780) [2ll3]. Later they confirmed the S+^(3520) state 
in two different decay modes (S+ —> KcK~'k^-, —)• pD^K~) at a mass of 3518.7 ± 1.7 MeV with an average 

lifetime less than 33 /s [3]. However, these states have not been observed either by BABAR [3, Belle 0111 in e“^e 
annihilation experiments or by LHCb at baryon-baryon collider experiments at CERN 0. In SELEX, production 
of doubly charm baryons with a large cross section occurred through baryon-baryon interactions which is totally 
inconsistent with fragmentation production. Further, the helicity angular distribution analysis also suggests that the 
pair of states, {S+(3443), S++(3460)} and {S+(3520), S++(3541)} form isospin doublets with isospin splittings 17 
and 21 MeV respectively. There is no precise understanding of these unusually large isospin splittings observed in 
the doubly charm baryons, unlike the small isospin splittings in the case of light and singly charmed baryons. An 
explanation may be possible if the coulombic electro-magnetic effect is much larger than strong interaction effect, 
leading to these baryons having a very compact size 0. In summary, the experimental status of doubly-charmed 
baryons has not been settled, however it is expected that consolidated analysis of large data collected from the ongoing 
experiments at LHCb and future experiments like PANDA ©FAIR, and Super Belle will shed further light. 

Doubly-charmed baryons are interesting systems as they provide a unique insight into the nature of the strong force 
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in the combined presence of both slowly moving heavy quarks along with the relativistic motion of a light quark. The 
excited spectra of these states and the splittings between them can help to understand how the collective degrees of 
freedom give rise to excitations in these systems. A comparison of these excitations with the corresponding spectra 
of singly and triply-charmed baryons, where the number of charm quark is one less and one more respectively, will be 
helpful to get information about quark-quark interactions. Doubly-charmed baryons are characterized by two widely 
separated scales: the low momentum scale, of order Aqcd, of the light quark and the relatively heavy charm quark 
mass. A doubly-heavy baryon can be treated as a bound state of a heavy antiquark and a light quark in the limit 
when the typical momentum transfer between the two heavy quarks is larger than Aqcd [101 111] . In this limit of 
quark-diquark symmetry, QQq o Qq, one can get definite prediction of spin dependent energy splittings [12[ig. It 
was argued mids] that because of heavy-quark symmetry, the doubly heavy baryons can be viewed as ultraheavy 
mesons, [QQ]q ^ Q'q and the hyperfine mass splittings in these systems are suppressed. Moreover, the doubly heavy 
baryon ground state of the form [QQ]j=iq will consist of chiral multiplets containing spin (1/2+, 3/2+) heavy spin 
fields [Hj. It is thus interesting to study these spin splittings to determine if the charm quark is sufficiently heavy 
to respect this quark-diquark symmetry. Doubly-charmed baryons have been studied over the years using various 
theoretical methods such as the non-relativistic [TMTB] as well as relativistic muD] quark models, heavy-quark 
effective theory mi, QCD sum rules |22{2S]i Feynman-Hellmann theorem mass formula |28j and Skyrmion 
model [22 • 

In light of existing and future experimental efforts to observe doubly-charmed baryons, it is desirable to have first 
principle predictions from lattice QCD. A quantitative description of the spectra of doubly-charmed baryons from 
the non-perturbative method of lattice QCD is valuable as it will enable a comparison between the lattice-computed 
spectra of doubly-charmed baryons to those obtained from potential models which have been successful for charmonia. 
Moreover, all results from such a first principles calculation will be predictions and thus naturally can provide crucial 
inputs to current and future experimental discovery. Given this significance of doubly-charmed baryons, it is desirable 
to study these states comprehensively using lattice QCD. Lattice QCD groups have studied the ground states of 
the doubly-charmed baryons using quenched NRQCD [30], quenched QCD with relativistic quarks [3TH33], and full 
QCD [5iMT] calculations. However, all previous lattice calculations involve only the spin 1/2 and spin 3/2 ground 
state spectrum of Sec and flee- It is expected that much more information about the interactions between two charm 
quarks and between charm and light quarks can be obtained by computing the excited state spectra of these baryons, 
including in particular the spin-dependent energy splittings, as well as by studying similar spectra for other spin-parity 
channels. Towards this goal, we report here the first attempt to compute the excited state spectra of doubly-charmed 
baryons using dynamical lattice QCD. The ground states for each spin-parity channel as well as their excited states 
up to spin 7/2 are computed and a few spin-dependent energy splittings are also studied. Similar studies of triply- 
charmed baryons have already been reported in Ref. |42| . and in a subsequent publication we will report results on 
singly charmed baryons. 

To extract the excited states of doubly-charmed baryons we follow the same procedure used in previous calculations 
for mesons [43ti42 and baryons [4^ I49H5T] . The gauge configurations we use are generated with 2-1-1 flavour clover 
fermions on anisotropic lattices [52E3]' For the charm quarks we also use clover fermions which are 0{a) improved 
at tree level in tadpole-improved perturbation theory. A large set of baryon operators which are first constructed 
in the continuum and then subduced into various lattice irreps |49| are used for this calculation. These operators 
transform as irreducible representations of SU{'A)f symmetry for flavor, S't/(4) symmetry for Dirac spins of quarks 
and 0(3) symmetry for orbital angular momenta. Baryon correlation functions are generated by using the “distillation 
method” |S1] and then variational fitting method is utilized to extract excited energies as well as to reliably determine 
the spins of these states. 

The layout of the paper is as follows. In the next section, we briefly describe lattice details, operator construction, the 
construction and analysis of correlation functions and our procedure for identifying the continuum spins of extracted 
states, which are already detailed earlier in Refs. IMMIIMI- In section III we present our results giving details of 
energy splittings in subsection IIIA. Finally, a summary of the work is presented in section IV. 


II. COMPUTATIONAL METHODS 

Over the last several years the Hadron Spectrum Collaboration (HSC) has adopted a dynamical anisotropic lattice 
formulation to extract highly excited hadron spectra. In this approach one uses a much finer temporal lattice spacing 
than in the spatial directions and exploits this higher resolution to extract highly excited states which decay rapidly at 
large Euclidean time separations with increasing noise to signal ratio. On the other hand, this avoids the computational 
cost that would come if one reduces the spacing in all directions. 


3 


A. The lattice action 

We used the tree-level Symanzik-improved gauge action and the anisotropic Shekholeslami-Wohlert fermion action 
with tree-level tadpole improvement and three-dimensional stout-link smearing of gauge fields for this work. The 
details of the formulation of actions as well as the method used to tune the anisotropic parameters can be found in 
Refs. [SH [S3]. In Table |l] we show the lattice action parameters of the gauge-field ensembles used in this work. We 


Lattice size 

atrrii 

atnis 

^cfgs 

m^/MeY mKlm^ atma 

16® X 128 

-0.0840 

-0.0743 

96 

391 1.39 0.2951(22) 


TABLE I: Details of the gauge-field ensembles used. Ndga is the number of gauge-field configurations. 

used the fl-baryon mass to determine the lattice spacing and obtained = 5.67(4) GeV. With an anisotropy of 
close to 3.5, this leads to Og = 0.12 fm, and total volume V ^ (1.9/m)^. We assume this volume is enough for a first 
study of doubly-charmed baryons. 

The charm quark action used here is similar to the light quark sector and the details are given in Ref. |47j . The rj^ 
meson ground-state was used to determine the bare charm quark mass. By studying the dispersion relation at low 
momenta, the action is made relativistic. As mentioned in 113, it is expected that the effects due to the absence of 
dynamical charm quark fields in this calculation will be small. 


B. Baryon operators 


Following Ref. |49j we construct a large set of baryon operators for doubly-charmed baryons. In summary the 
construction has two steps: a set of continuum operators with well-defined continuum spin is found and then these 
operators are subduced to the irreducible representations (irreps) of the octahedral group on the lattice. A set of 
continuum baryon interpolating operators with well-defined continuum spin are constructed as 


1 ~ (g) S-^s 0 


( 1 ) 


where J-, S and V represent flavor, Dirac spin and spatial structure respectively while the subscripts Si specify 
the permutation symmetry in the respective subspaces. The details of these permutation symmetries and their 
combinations are given in Refs. |49l 150] . 

By permutation symmetry, one can argue that the flavor structure of QQq and QQs will be the same as qqs 
combinations, which are the flavor structures of light S baryons. Hence the possible flavor-symmetry structures for 
doubly charm baryons are totally symmetric (S) in all three flavor labels belonging to the decuplet flavor constructions 
(IOf) as well as symmetric (MS) and antisymmetric (MA) in the first two flavor labels belonging to the octet flavor 
constructions (8f)- For the decuplet structure, the flavor labels being symmetric, the remaining spin and spatial part 
should be combined symmetrically to form an overall symmetric interpolating operator. For the octet structure, the 
flavor structure being MS and MA, the symmetry in the remaining part of the baryon operator, excluding the color 
labels, should be MS and MA to form a symmetric interpolating operator. 

The spatial- as well as spin-symmetry combinations used are also presented in Ref. [49j . Up to two covariant 
derivatives are considered and combined so as to transform as orbital angular momentum, L, with maximum accessible 
values of 0, 1 and 2. The subset of operators formed by considering only the upper two-component of the four- 
component Dirac-spinor in the Dirac-Pauli representation are labelled non-relativistic, while all other operators are 
called relativistic. Another subset of operators with D = 2 and L = 1 in the mixed symmetric and mixed antisymmetric 
combinations are identified as hybrid operators because of their essential gluonic content; these operators vanish in the 
absence of a gluon field. In Table jlj we show the allowed spin-parity patterns based on non-relativistic quark spinors 
for up to two covariant derivatives. Note that with the non-relativistic operators, it is not possible to construct 
a negative-parity state beyond spin | , even with non-local operators using two derivatives. Use of relativistic 
operators along with non-relativistic ones enable us to extract higher negative-parity states as well as higher excited 
states. These continuum operators are then subduced to the irreps of the cubic group. The three irreps of the double¬ 
valued representations of the octahedral group for half-integer spins are Gi, G 2 and H. The details of this subduction 
procedure to obtain the lattice operators was discussed in Ref. |49|. In Table III we show the number of operators 


that we obtained after subduction and are used for this study. Both the number of positive (g), negative (u) parity 
operators as well as the number of operators with non-relativistic quark spinors (NR) and with hybrid content are 
shown in this table. 
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TABLE II: Allowed spin-parity patterns based on non-relativistic quark spinors. For each covariant derivative (D = 0, 1 and 
2), the SU(3) flavor, the quark spin S and orbital angular momentum L are listed followed by the allowed values. The total 
number of operators is listed as Ni for derivative i. The lower part of the table gives the two-derivative hybrid operators based 
on non-relativistic quark spinors, as discussed in m- 



Gi H G2 

g u g u g u 

Total 

55 55 90 90 35 35 

Hybrid 

12 12 16 16 4 4 

NR 

11 3 19 4 8 1 


TABLE III: The total number of operators, with covariant derivatives up to two, obtained after subduction to various irreps. 
The number of non-relativistic (NR) and hybrid operators for each irreps and for both parities are also mentioned. 


C. Analysis of baryon correlators : Distillation and variational methods 

After constructing a large set of operators for each irrep A, we calculate the matrix of correlation functions 

C,^{t = tf- U) = {0\OSf)O]iU)\0). (2) 

between a baryon source at time ti and sink at time tf. The distillation method |54j provides an efficient means 
of constructing the large correlation matrices needed for this analysis. For this work, the method was realized by 
constructing the distillation operator with 64 eigenvectors of the gauge-covariant Laplacian and then correlation 
functions are computed from 4 time sources. The matrix of correlation functions was then analyzed using the 
variational method, which proceeds by solving a generalized eigenvalue problem of the form 


( 3 ) 
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here C(t) is the matrix of correlators at time-slice t, i.e. Cij{t) = {Oi{t)Oj{0)), A„(t,to) are the principal correlators 
and n"’s are the eigenvectors which determine the overlap of an operator to a particular state. The details of our 
fitting procedure is given in Ref. [33 and was utilized to extract the excited state spectra of mesons as well 

as of baryons [31I3SHSI]- In Figure [2 and Figure]^ we plot some examples of fits to the principal correlators for flee 
baryons in the Hg irrep and for Sec baryon in the G 2 g irrep. 



FIG. 1: Principal correlator fits for six states of flee baryons in irrep. Hg that are identified as J = 3/2+. Data points are 
obtained from Fits are carried out using a fitting form X„{t) = (1 — *■*“*“! with three 

fit parameters, m„ : the state that we are looking for, m/ : where all other excited states are grouped together, and An : 
which is related to overlap factor. The lines show the fits and one sigma-deviation according with to = 12; the grey points are 
not included in the fits. 


D. Rotational symmetry 

We constructed operators in the continuum which were then subduced onto the lattice irreps to form lattice 
operators. This allows us to determine the properties of these states in the continuum with some confidence. Following 
Refs. [42l I4M5T] . Figured shows the normalized correlation functions, Cijj^CuCjj, for 5cc baryons transforming 
under Hg and determineefat time-slice 5. The normalization ensure diagonal entries are unity and off-diagonal entries 
are less than 1. Various operators are represented by following abbreviations : non-relativistic (n), relativistic (r), 
non-hybrid (I) and hybrid (2). There are 90 operators used in this irrep, including operators up to two derivatives. 
The solid lines divide these operators into spins |, | and and the dashed lines separates operators defined above. 
As is evident from the figure, the matrix is close to block diagonal implying that there are small correlations between 
operators subduced from different continuum spins. This suggests that there remains a remarkable degree of rotational 
symmetry in the matrices of correlation functions obtained for these lattice operators, similar to the light, strange 
and triply-charmed baryons [421150| . Similar block-diagonal matrices of correlation functions are observed in other 
irreps and for flee baryons as well. 

From the correlators one can also identify the flavour mixing between different operators. For example, it is evident 
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FIG. 2: Same as Figure]^ but for baryons in irrep. G 2 g- 


from Figure that there is strong mixing between relativistic operators belonging to octet and decuplet. However, 
this flavour mixing is less evident for non-relativistic operators. As was observed in our previous study on triply- 
charmed baryons |42j . we found that there are additional suppression in mixing for non-relativistic operators with a 
given J, but with different L and S, compared to those with the same J, as well as the same L and S. However, 
this suppression is not present for relativistic operators. On the contrary, in most of the cases involving relativistic 
operators, we found that mixings are comparable for operators with given J, but with different L and S to those of 
operators with the same J, as well as same L and S. 

1. Continuum spin identification 

One main difficulty of lattice calculations of spectra is the identification of the spin of an extracted state at finite 
lattice spacing, particularly for the state onto which various operators from different irreps contribute. In Ref. |43j a 
well defined procedure is adopted for spin identification by using the overlap factor. The overlap-factors of an operator 
Oi to a state n are Zf = (n|O||0) which can be shown [43] to occur in the spectral decomposition of the matrices of 
the correlation functions as, 

^n* 2^n 

corra,{t)=Y,^^e-^-K (4) 

„ 2m„ 

One can use the orthogonality for the eigenvectors u”^C'(to)'c™ = J"’™ to show that the overlap factors can be obtained 
from the eigenvectors using the relation 

Zf = (5) 

In Figure 4 we show an array of histograms of the normalized overlap factors, Z = ^ operators 
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FIG. 3: The normalized correlation matrix, Cij f ^/CiiCjj, at t/at = 5 are shown for Sec baryons in the Hg irrep., according 
to the darkness scale at the side. The operators are ordered such that those subduced from spin 3/2 appear first followed by 
spin 5/2 and then spin 7/2. The correlation matrix is observed to be mostly block diagonal in terms of spins which signifies 
the rotational symmetry on our lattice. There is significant mixing between octet (8f) and decuplet (IOf) operators which is 
not present in the corresponding plot for light quark E(uus) having same operator structures (see figure 2 of Ref. [50]). This 
signifies that the SU{3) flavor symmetry is badly broken for doubly-charmed Sec baryons. Non-relativistic (n), relativistic (r), 
non-hybrid (1) and hybrid (2) operators are also identified for each of those spins. The strength of flavor mixing is maximum 
between relativistic operators (r), while for non-relativistic operators this mixing is relatively mild. 


onto some of the spin identified lower-lying states in each of the lattice irreps. The normalization Z = \z’'] 

such that the largest overlap factor of that operator across all states is unity. 

As in Ref.jSD], a “matrix” plot of the overlap factors can be used to depict the dominant contribution to the 
low-lying states from each operator. As examples, Figure and Figure show such plots of the normalized overlap 
factors, Z”, of an operator i to a given state n, as defined by eq. ([^, for Sec and flee respectively. All values of 
are normalized as above. Figure shows that state 0, the ground state, and excited states 1, 2, 3, 4, 9, 12 and 15 are 
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FIG. 4: Histogram plot showing normalized overlap factor, Z, of a few operators onto some of the lower-lying states in each 
of the lattice irreps. Z’s are normalized according to , so that the largest value for that operator across all states is 

equal to unity. Top row is for irrep Gig, middle one is for Hg, and the bottom one is for G 2 g irrep. Black bars correspond to 
spin-1/2 operators, red for spin-3/2, green for spin-5/2 and blue represents spin-7/2 operators. Lighter and darker shades on 
the top of every bar represent the one sigma statistical uncertainty. 


naturally identified as = I"*" states while states 5, 7, 8, 10 and 13 are identified as = I"*" states and similarly, 

states 11 and 14 can be identified as = I"*" states. However, in order to confirm the reliability of the identification 
of a state with a given spin greater than 3/2 we compare the magnitudes of overlap factors of a particular operator 
from different irreps. This will be discussed later. These plots help us to identify the structure of a state from the 
types of operators which construct it. For example, in Figure]^ the spin-3/2-|- ground state has predominant overlaps 
from the non-relativistic non-hybrid as well as relativistic non-hybrid type operators. Similarly, the states 11 and 14, 
which are identified as I”*" states, have predominant overlaps to non-relativistic non-hybrid operators. Strong hybrid 
content was observed for a number of states by identifying their strong overlaps with hybrid interpolating operators. 

To identify the spin parity of a state we followed the same method detailed in [43] and used for light mesons [441146] , 
baryons gniisn], charm mesons m as well as heavy-light mesons [JS] calculations. To identify a spin that subduces 
into a single irrep is relatively straightforward and studying the overlap factors used in histogram and matrix plots 
we can identify spin-^ and spin-| states. For spin-| and spin-| states, which are subduced into multiple irreps 
overlap factors from different irreps must be compared. As the continuum limit is approached, the factors of a given 
continuum operator to a particular state obtained from various subduced irreps should be the same. For a fine lattice 
spacing these overlap factors should thus be close to each other. For example, the spin-7/2 continuum operator, 
(3/2+)i_s 0 subduced to irrep Gig, Hg as well as to G 2 g- This near degeneracy of overlap factors can 

be used to identify this state as spin-|^ state. In Figure [^we compare a selection of Z-values for states conjectured 

to be J = 1^ (top two plots), I"*" (middle two plots), | (bottom left) and | (bottom right) which appear over 
multiple lattice irreps. The continuum operators considered are noted along the lower edge of each plot. Z-values 
obtained for a given operator but from different irreps are found to be consistent with each other which helps us to 
identify the spin of these given states. 

After identifying the spin of a state with matching overlap factors we check whether the energy of this state 
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FIG. 5: “Matrix” plot of the normalized overlap factor, Z", of an operator i to a given state n, as defined by eq. Z'^ 
are normalized according to ^ax\zv-] ’ that for a given operator the largest overlap across all states is nnity. This plot 
corresponds to the extracted states of Hcc, in the Hg irrep. Darker pixels indicate larger values of the operator overlaps as in 
Figure]^ Various type of operators, for example, non-relativistic (n) and relativistic (r) operators, as well as non-hybrid (1) 
and hybrid (2) operators are indicated by column labels. In addition, the continuum spins of the operators are shown by 3/2, 
5/2 and 7/2. State 0, the ground state, and excited states 1, 2, 3, 4, 9, 12 and 15 could be identified as = I”” states by the 
overlap to various types operators according to pixel strengths. States 5, 7, 8, 10 and 13 could be identified as states; 

similarly, states 11 and 14 could also be identified as states. 


determined over different irreps also matches. To achieve this check, the corresponding principal correlators across 
the irreps are simultaneously fit to a single energy. In Figure these joint fits are shown for the principal correlators 
obtained from three different irreps for a spin-7/2+ state (top plot) and from two different irreps for a spin-5/2+ state 
(bottom plot). 


III. RESULTS 


In this section, our results for the doubly-charmed baryons spectra with spins up to | and with both parities are 
presented in terms of energy splittings. In general, this tends to reduce systematic uncertainties in lattice calculations, 
including those from the scale setting procedure used. For this work, the charm quark mass parameter in the lattice 
action was determined by ensuring the physical value and lattice estimate of the rjc meson mass |47j agreed, once the 
lattice spacing is determined using the U-baryon. The 77 c meson has the same number of charm quarks as the doubly- 
charmed baryons. Hence, we show the spectra of Sec and Hec baryons with the mass of the rye meson subtracted. 
These energy splittings expose the binding energies of the extra light and strange quarks inside these baryons. A few 
energy splittings such as the hyperfine splittings between the extracted states are also computed. 

As mentioned earlier, the experimental status of doubly heavy baryon discovery is uncertain. It is thus important 
to compare the ground state spectra of the two doubly heavy baryons obtained from different calculations. In Figure 
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FIG. 6: Same as Figure]^ but for Slcc baryons in the irrep Hg. 


the only experimental (SELEX) result and other lattice as well as various model results. In Figure 9 we show similar 
results for flee baryons. Our results are at pion mass 391 MeV, results for ETMC [35], PACS-CS |^T , Bali et. al [30] . 
and Briceno et. al m are extrapolated to the physical pion mass, while ILGTI |33| results are at pion mass 390 
MeV. While the lattice spacing in the temporal direction (at) for this work is 0.0351 fm, for ETMC, at = 0.056 fm, 
for PACS-CS, at = 0.0899 fm, for Bali et. al, at = 0.0795 fm and for ILGTI, at = 0.0582 fm. Results for Briceno et. 
al are extra pola ted to the c ont inuum limit. 

In Figure 11 and Figure 12 we show the spin identified full spectra, up to J = | extracted from our lattices of 


Sec and rice baryons, respectively. The states inside the magenta boxes are those with relatively large overlap onto 
non-relativistic operators and the states with thick borders corresponds to the states with strong hybrid content as 
defined in Ref. m- 

Note that in the lowest two positive parity bands and the lowest negative parity band, the number of states for 
each spin agrees with the expectation shown in Table [llj That table gives the number of allowed quantum numbers 
by SU{6) X 0(3) symmetry for operators with up to two derivatives (D). For and = I"*", the number of 

allowed quantum numbers is 7 (1 from D = 0 and 6 from D = 2) and 9 (1 from D = 0 and 8 from D = 2) respectively. 

For other positive parity quantum numbers with = 1^ ^-'^d these numbers are 5 and 2 respectively. The left 
sides of Figure 11 and Figure 12 show positive parity states. The numbers of states in the lowest two bands (first 
two and inside box) match exactly with the allowed quantum numbers mentioned above. Note that while we use 
the full set including both non-relativistic and relativistic operators not in Table [^ we still obtain the same number 
of states allowed with only non-relativistic operators. Similarly, for negative parity with ^ , | , | , the 

allowed number of quantum numbers are 3, 3, 1 and 0, respectively. On the right sides of Figure 11 and Figure [T^ 
the lowest band (inside box) also has exactly the same number of states. This agreement of the number of low lying 
states between the lattice spectra obtained in this work and the expectations based on non-relativistic quark spins 
implies a clear signature of SU{6) x 0(3) symmetry in the spectra. Such S17(6) x 0(3) symmetric nature of spectra 
was also observed in Ref. [42l [50]. Note there are no negative parity spin-7/2 state in that table and the negative 
parity spin-7/2 state is obtained from the inclusion of relativistic operators. We are also able to identify one state 
with strong overlap onto a hybrid operator. Though there are more quantum numbers accessible with the operators 
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and I (bottom right). The operators used are mentioned at the bottom of each plot. Z-values obtained for a given operator, 
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FIG. 8: Joint fit of the three principal correlators : top figures - for a representative spin-|^ state from three different irreps, 
and bottom figures - for a spin-|^ state from two different irreps. Joint fits with X^/dof ~ 1.9 (for both cases) provide the 
masses of these states. 


in Table |n] with hybrid structures, we could not clearly identify those. This is because, as noted in Ref. m , it is not 
meaningful to interpret the higher excited states in terms of SU{6) x 0(3) symmetry. 

Figure shows the extracted spectra with a different subtraction. The doubly-charmed baryons have two 
charm quarks and one light quark. On the other hand, D meson has one charm quark and one light quark. The 
computed energy splittings between Sec and the ground state of D meson are shown in Figure In the same figure 
the energy splittings between Vice and the ground state of Dg are shown. For each spin, the left column indicates 
splittings of Sec — D while the right column shows Vice — Dg splittings. The inset at the top hgures shows the positive 
parity ground states of spin-1/2 and spin-3/2 baryons. Results for D and Dg mesons are taken from Ref. [JB]. By 
subtracting the D and Dg meson masses from Sec and Vice baryons we effectively leave only the energy of the excitation 
of a single charm quark in both cases. One would thus naively expect that both spectra will be equivalent. This is 
almost true for the lowest state in each spin parity channel, except for negative parity spin-7/2 state, as shown in 
Figure [T^ That state is obtained from relativistic operators and it is expected that such naive expectation may not 
hold there. This is also seen for the excited states where the contribution from relativistic operators is much larger. 

As in our previous studies in charmonia [47], charmed-strange mesons [48] and triply-charmed baryons [42], the 
systematic uncertainty due to 0(a) discretisation artefacts was investigated in another calculation in which the spatial 
clover co-efficient was boosted from the tree-level Cg = 1.35 to Cg = 2. At this value our extracted hyperfine splitting 
was found to be physical [47]. The mass splittings were found to differ by around 40-45 MeV between the two 
calculations. In this study also we found a similar mass difference between two calculations, indicating an equivalent 
scale for the uncertainties in this calculation. 
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FIG. 9: Ground state masses of spin-1/2 and spin-3/2 doubly-charmed 11 baryons as a splitting from rjc meson mass. Our 
results are shown by the red filled circle (HSG). Other lattice as well as model results are also shown. 


A. Energy Splittings 


The energy splittings between various excitations in a spectrum provide important information about the nature of 
interactions needed to excite those states. The energy splittings also provide inputs for building models to describe 
these states successfully. In Figure we show energy splittings of the ground states of each spin parity channel from 
the lowest state in that parity channel. For the positive parity, the lowest state is and for negative parity 

the lowest state is ^ . It is interesting to note that both for S^c and flee this splittings are almost same. This 

indicates that the interquark interactions, which are responsible for these splittings, are similar in these two different 
type of baryons. 

The most notable spin dependent baryon energy splitting to consider is the hyperfine splittings between | and ^ 
states, for example splitting between A and nucleon. For doubly-charmed baryons we also compute this splitting and 
show in Figure 15 for Sec (left plot) and flee (right plot) baryons. Our results (red circles) are compared with other 
lattice results (Hue squares) as well as with various model results. Note our results for Sec are for pion mass 391 
MeV. Results for ETMC [38], PACS-CS [41], Bali et. al [40], and Briceno et. al [37] are extrapolated to the physical 
pion mass, while ILGTI [3S] results are at pion mass 390 MeV. 

In Figure [^ we compare these hyperfine splittings for Sec and Occ baryons as a function of quark mass. In the 
ai-axis of that figure we use the square of the pseudoscalar meson mass {rrips) while y-axis shows hyperfine splittings 
at those pseudoscalar meson masses. Along with Sec and flee we also show splittings between spin-| and spin-^ states 
of riccc baryon. The later spin dependent splitting is not hyperhne in nature as both of them are of decuplet type. As 
in Ref. [22], for the positive parity flccc baryons we take the spin-orbit splitting between i? 3 (|''') and states 

which have same L and S values. For the negative parity flccc baryons we take spin-orbit splitting between i?o(| ) 
and Ao(| ) states which also have same L and S values. 


In Ref. 


we evaluated energy splittings between triply-flavoured baryons^ and studied their quark mass depen- 


by triply-flavoured baryons, we mean A , Qcce and baryons. 
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dence. It was observed that various splittings decrease significantly with quark masses. Following heavy quark effective 
theory one can expand the mass of a heavy hadron, with n heavy quarks, as = tiMq+A+B/ mQ+0{l/mQ2 ) [55] . 

With that expectation, we fitted various energy splittings of triply-flavoured baryons with a form a + b/rups, with nips 
the pseudoscalar meson mass, and obtained good fits. Motivated by that, here we also fit similar energy splittings for 
doubly-charmed baryons. There is a however a difference for doubly-charmed baryons because to show quark mass 
dependence one needs to use triply-charmed baryons data also. For the triply-flavoured baryons in Ref. [42], states 
at different quark masses are obtained from the same operators and so are simple to compare. However, in this work, 
doubly- and triply-charmed baryons are not obtained from the same operators except for positive parity spin 3/2 and 
spin-7/2 states. We thus consider energy splittings only for these states. More specifically, we calculate the following 
energy splittings : El^{ccu)—Du{cu) , ^l*^{ccs) — Ds{cs) and Q.*^^{ccc)—r]c{cc), and E*^{ccu)—D^{cu) , H*^(ccs)—i4*(cs) 


and H*^^(ccc) — J/il){cc) and plot in Figure 17 As in Ref. |32|, we fit these splittings to the form a + b/mps and obtain 


reliable agreement. Energy splittings from vector mesons are also fitted with a constant. We emphasis that energy 
splittings at very light quark masses should not be calculated using the above form, which is valid only for heavy 
quarks. 

It is interesting to note that one can extrapolate the fitted results to the bottom mass to obtain the energy splittings 
of fl/cb(§''') ~ and ~ at the mass of meson. The difference between these two quantities yields 

the mass splitting of B* — Be- Using this, we obtain the mass splitting of B* — Be as 80 ± 8 MeV (fitting with a 
form a + b/rups) and 76 ± 7 MeV (fitting with a constant term). This result agrees very well with those obtained by 
potential models [56] . It does not however agree with other lattice QCD prediction [571 - 159] which provide much lower 
values. It will be interesting to study the effects of the lattice cut-off on the hyperfine splittings of these hadrons 
with the possibility of future discovery of B* meson. Only one charmed-bottom meson has been discovered so far [T] . 
Assuming it as a pseudoscalar meson with mass 6277 MeV and using the extrapolated value of ~ we 

empirically predict the mass of U/^j,(3/2+) to be 8050 ± 10 MeV which is consistent with the results from various 
models cniiiiiiiiiii and a recent lattice calculation [? ]. 
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FIG. 10: Ground state masses of spm-1/2 and spin-3/2 doubly-charmed H^c baryons as a splitting from r]c meson mass. Our 
results are shown by the red filled circle (HSG). Our u quark mass is not physical and at pion mass 392 MeV, results for ETMC 
[38], PAGS-GS [41], Bali et. al [40], and Briceno et. al [37) are extrapolated to the physical pion mass, while ILGTI |39) results 
are at pion mass 390 MeV and 340 MeV respectively. 
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FIG. 11: Spin identified Hcc baryon spectra for both parities and with spin up to 7/2. Energy splittings of the H^c baryons 
from the mass of m,,^ meson, which has same number of charm quarks, are shown here. The states inside the pink boxes are 
those with relatively larger overlap to non-relativistic operators and the states with thick borders corresponds to the states 
with strong hybrid content. The number of states inside these boxes matches with the expectations based on non-relativistic 
quark spins, as shown in Table |IT] This agreement of the number of low lying states between the lattice spectra obtained in 
this work and the expectations based on non-relativistic quark spins implies a clear signature of SU{6) x 0(3) symmetry in the 
spectra. 



FIG. 12: Spin identified spectra of Qcc baryons for both parities and with spin up to 7/2. Energy splittings of the Qcc states 
from the mass of r/c meson, which has same number of charm quark, are shown here. All other details are same as in Figure 
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FIG. 13: Energy splittings of H^c and lice baryons from D and Ds mesons, respectively. Splittings for the positive parity 
baryons are in the left plot and those for negative parity baryons are shown in the right plot. For each spin, left column is for 
splittings of Hcc(ccm) baryons from the ground state D{cu) meson and the right column is for splittings of flcc(ccs) baryons 
from the ground state Ds{cs) meson. The inset at the top figures are for the positive parity ground states of spin-1/2 and 
spin-3/2 baryons. Results for D and Ds mesons are taken from Ref. |48 |. 
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FIG. 14: Energy splittings (in units of GeV) of the ground states of each spin parity channel from the lowest state in that 
parity channel. Splittings for Hec are shown with red color rectangle and those of are represented by blue color shredded 
rectangle. For the positive parity, the lowest state has and for the negative parity the lowest state has | . 


IV. CONCLUSIONS 

In this work, results from the first non-perturbative calculation on the excited state spectroscopy of the doubly- 
charmed baryons with spin up to 7/2 are presented. We performed our calculations using lattice QCD with dynamical 
2-1-1 flavours clover quarks on anisotropic lattices. There is no clear experimental observation to date of a doubly- 
charmed baryon. It is important for theory to determine if various model results using quark-diquark symmetry or 
otherwise, can be tested with a non-perturbative method for these systems for which the interaction is governed by 
two widely separated scales. Our extensive study of doubly charmed baryons has made a first attempt to connect 
possible experimental and theoretical studies. 

We employ a large basis of creation operators and use the distillation technique to calculate the doubly-charmed 
baryon correlation functions. A variational fitting method was employed to extract the spectrum. We observe 
approximate rotational symmetry for these operators at the scale of hadrons. Having realized a good rotational 
symmetry we determine the overlap factors, as in previous studies [^^51) . We are able to extract states reliably with 
spin up to 7/2 and also studied the mixing between various operators. 
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at pion mass 392 MeV, results for ETMC [38], PACS-CS [41] . Bali et. al [40], and Briceno et. al [32] are extrapolated to the 
physical pion mass, while ILGTI |39| results are at pion mass 390 MeV and 340 MeV respectively. 
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FIG. 16: Mass splittings between spin-|^ and spin-|^ states (left figure) and between spin-| and spin-| ground states 
(right figure) are compared for Ecdccu), ^lcc{ccs) and flccc(ccc) baryons. For Uccc, mass splitting is between iJo(l/2''") and 
753 ( 3 / 2 ^) states, which is actually due to spin-orbit coupling, while for Sec and flee, these are between respective ground states 
and are due to hyperhne splittings. 


The main results are shown in Figure E] and Figure [T^ for ^cc and flee baryons, respectively. The ground states 
for spin 1/2 and 3/2, which have more relevance to recent experiments, are shown in Figure and Figure 10 As 
in Ref. g5|[S0], we also find bands of states with alternating parities and increasing energies. We also observed the 
number of extracted states of each spin in the three lowest-energy bands and the number of quantum numbers expected 
based on weakly broken SU{6) x 0(3) symmetry agree perfectly, i.e., the doubly-charmed baryon spectra remarkably 
resemble the expectations of quantum numbers from non-relativistic quark model |63H65j . This symmetry was also 
observed for light, strange m and triply-charmed baryon spectra [12]. For positive parity states this agreement does 
not get spoiled even with the inclusion of non-relativistic hybrid operators. However, it is expected that this band 
structure will persist with the inclusion of relativistic operators which contribute more to the higher excited states. 
The extracted spectra for the higher excited states also support this observation. Note our operator set does not 
include any multi-hadron operators. It is expected that inclusion of those operators, particularly those involving light 
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FIG. 17: Energy splittings between positive parity spin-3/2 baryons and psendoscalar meson as well as vector mesons are 
plotted against the square of the pseudoscalar masses. Left figure is the splittings of : H*c(ccu) — (cu), fi*^(ccs) — Ds (cs) and 
flccc(ccc) — rjc{cc), and right figure includes following splittings : S*c(ccu) — D^(cu), Q*c(ccs) — Dl{cs) and Q*^^(ccc) — Jlip(cc); 
they are plotted against the square of the pseudoscalar masses (i.e., at Du , Ds) and rjc mass. We fit the quark mass dependence 
with a form a + hjm^s (right figure is also fitted with a constant term). The fitted results are shown by solid lines with shaded 
regions as one sigma errorbars. 


quarks, may modify some of these conclusions. We are also able to decode the structure of operators leading to a 
particular state: whether constructed by relativistic, non-relativistic, hybrids, non-hybrid types or a mixture of them 
all. However, this identification is not possible for negative parity states and highly excited positive parity states, as 
argued in Ref. [50] . 

Our calculated spectra do not support the chiral multiplet structure of doubly heavy baryons speculated in Refs. 
[HIS]. The ground states of doubly charmed baryons are not degenerate and both for S^c and flee hyperhne splittings 
are around 80 — 100 MeV, as shown in fig 15. The calculated spectra do not match the diquark picture and instead 
they are remarkably similar to the expectations from non-relativistic quark models with an SU(6)x 0(3) symmetry. 

The study of the energy splittings between various excited states is quite helpful in revealing the nature of interquark 
interactions and a detailed knowledge of them could help to build successful models. We calculated hyperhne mass 
splittings between spin-| and spin-^ states of .rice and flee and compared those with various other lattice and 
model predictions. 

Energy splittings between the ground states of different spins were also evaluated and we observed that the hierarchy 
of the hrst few energy excitations in and flee are quite similar. This indicates the involvement of similar dynamics 
to excite these states. To study the quark mass dependence of the energy splittings we compared results for Sec (ecu), 
f2cc(ccs) and Hccc(ccc) baryons for which there is a common ‘cc’ diquark and a varying quark from light to charm. 
Encouraged by a successful htting of the mass splittings in triple-havoured baryons |42j we studied similar mass 
splittings of doubly-charmed baryons. Here also we hnd that a heavy quark motivated form a + bjmps can ht quite 
successfully energy splittings like : Scc(ccu) — Du{cu) ,fll^{ccs) — Ds{cs) and fll^^{ccc) — ridcc), and 'E.*^{ccu) — 
Du{cu ), n*c(ccs) — Dl{cs) and fl*^^(ccc) — Erom the Htted results we are able to predict R/ — He = 80 ± 8 

MeV and n;ef,(3/2+) = 8050 ± 10 MeV. 
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